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Abstract 

Using the D0 detector at the 1.8 TeV pp Fermilab Tevatron collider, we 
have measured the inclusive dijet mass spectrum in the central pseudorapidity 
region |r?j c t| < 1-0 for dijet masses greater than 200 GeV/c 2 . We have also 
measured the ratio of spectra a(\rjj et \ < 0.5)/<r(0.5 < |r?j ct | < 1.0). The order 
a 3 QCD predictions are in good agreement with the data and we rule out 
models of quark compositeness with a contact interaction scale < 2.4 TeV at 
the 95% confidence level. 
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High transverse energy (Et) jet production at a center of mass energy of 1.8 TeV probes 
the structure of the proton down to a distance scale of 10~ 4 fm. A measurement of the dijet 
mass spectrum can be used to verify the predictions of quantum chromodynamics (QCD) 
for parton-parton scattering and to constrain the parton distribution functions (pdf) of the 
proton. Additionally, new physics such as quark compositeness |1[ would be revealed by an 
excess of events in the dijet mass (M) spectrum at high masses with respect to the predic- 
tions of QCD. A previous analysis by the CDF collaboration of the inclusive jet cross section 
reported an excess of jet production at high Et- More recent analyses of the dijet angular 
distribution by D0 || and CDF || have excluded at the 95% confidence level models of 
quark compositeness in which the contact interaction scale is below 2 TeV. Most recently, 
an analysis of the inclusive jet cross section by D0 || shows good agreement between the 
theory and data. This paper presents a new improved measurement by D0 of the inclusive 
dijet mass spectrum (uncertainty reduced by a factor of four at M = 200 GeV/c 2 relative to 
the previous CDF measurement ||) and improved limits on the contact interaction scale. 

The outgoing partons from the parton-parton scattering process hadronize to form jets 
of particles. These jets were identified in the D0 detector using uranium/liquid- argon 
calorimeters which cover a pseudorapidity range of 1 77 1 < 4.1 (rj = — ln[tan(0/2)], where 
6 is the polar angle relative to the proton beam direction). The calorimeters have a jet 
transverse energy resolution of 10% (5%) for E T of 50 (300) GeV. 

Events with at least one inelastic interaction during a beam crossing were identified using 
scintillator hodoscopes, and the primary event vertex was determined using tracks recon- 
structed in the central tracking system. Event selection occurred in two stages. First, a 
minimum transverse energy was required in a region (A77 x A0 = 0.8 x 0.6) of the calorime- 
ter. Jet candidates were then reconstructed online with a cone algorithm of opening angle 
1Z = 0.7 in rj-(f> space (<ft is the azimuthal angle), and the event was recorded if any jet 
Et exceeded a specified threshold. During the 1994-95 run, the thresholds were 30, 50, 
85, and 115 GeV, with integrated luminosities of 0.353±0.027, 4.69±0.37, 54.7±3.4, and 
91.9±5.6 pb _1 , respectively. The luminosities of the 30 and 50 GeV triggers were deter- 
mined by matching their dijet cross sections to that measured for the 85 GeV trigger. This 
resulted in an additional uncertainty of 4.9% in the luminosities of the 30 and 50 GeV trig- 
gers. 

Jets were reconstructed offline using an iterative jet cone algorithm with 1Z = 0.7 ||. 
Jet Et is defined as the sum of the Et in each cell within the cone. The jet was centered 
on the -&r-weighted pseudorapidity and azimuth of the jet. The jet Et and direction were 
then recalculated until the cone direction was stable. If two jets shared energy, they were 
combined or split based on the fraction of energy shared relative to the Et of the lower- Et 
jet. If the shared fraction exceeded 50%, the jets were combined and the direction recalcu- 
lated. Otherwise, the jets were split and the energy in each of the shared cells was assigned 
to the nearest jet. The directions of the split jets were then recalculated ||. 

A significant fraction of the data was taken at high instantaneous luminosity, which re- 
sulted in more than one pp interaction in a beam crossing leading to an ambiguity in selecting 
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the primary event vertex. After event reconstruction, the two vertices with the largest track 
multiplicity were retained. The quantity St —\T,E^\ was calculated for both vertices, and 
the vertex with the smaller St was selected. The uncertainty on the mass spectrum due to 
the choice of vertex was 2%. The vertex was required to be within 50 cm of the detector 
center. This cutoff was 90 ± 1% efficient, independent of the dijet mass. 

Backgrounds from noise, cosmic rays, and accelerator losses were reduced to an insignif- 
icant level by applying jet quality criteria. For an event to be accepted, the two leading- E T 
jets were required to satisfy these quality criteria. Contamination from backgrounds was 
< 2% based on Monte Carlo simulations and visual inspection of dijet events of high mass. 
The overall jet selection efficiency for 1 77 1 < 1.0 was measured as a function of jet Et, giving 
97 ± 1% at 250 GeV, and 94 ± 1% at 400 GeV. 

The transverse energy of each jet was corrected for the underlying event, additional inter- 
actions, calorimeter noise, the fraction of particle energy that showered hadronically outside 
of the cone, and for the hadronic response |J. At 77 = 0, the mean total jet energy correction 
was 16% (12%) at 100 GeV (400 GeV); the correction uncertainty was less than 2.5% of the 
jet E T . 

For each event that passed the criteria, the dijet mass, M, was calculated, assuming that 
the jets are massless, using M 2 = 2E^ tl E^ t2 [cosh.(Ai]) — cos(A0)]. 

The steeply falling dijet mass spectrum is distorted by jet energy resolution (and 
to a negligible extent by r\ resolution). The dijet mass resolution was calculated us- 



ing the measured single-jet resolutions and the pythia pi Monte Carlo event gener- 
ator. This resolution depends on the Et and rj distributions of the two leading Et 
jets in each event. The observed mass spectrum was corrected with an ansatz function 
F(M') = BM'~ a [l — (M' I \fs)]~P convoluted with the mass resolutions, to obtain the 
smeared ansatz f(M) = / °° F(M')p(M' - M, M')dM' (where p is the mass resolution) , such 
that the number of events in any given mass bin i is given by integrating / over that bin. 



The data were then fitted using a binned maximum likelihood method and the minuit JTT 
package to determine the values of B, a, and (3. The unsmearing correction for each mass 
bin is then given by C t = J FdM/ J fdM (Q = 0.96 (0.92) at 210 (900 GeV/c 2 )). 



The dijet mass spectrum was calculated using 

K = d 3 a/dMdr ]l dr] 2 = {N.d) / (C t eAM A^Arfr), 

where Ni is the number of events in mass bin i, Li is the integrated luminosity for that bin, 
e is the efficiency of the vertex selection and jet quality cuts, AM is the width of the mass 
bin, and A^ lj2 is the width of the rj bin for jets 1 and 2. The spectrum was calculated for the 
pseudorapidity range | i]- ]et \ < 1.0 (where both jets are required to satisfy the 77 requirement), 
in mass ranges starting at 200, 270, 350, and 550 GeV/c 2 , corresponding to the software jet 
E T thresholds of 30, 50, 85, and 115 GeV. 
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The cross section for the mass spectrum is plotted in Fig. [I], and given in Table |. 
The data are plotted at the mass-weighted average of the fit function for each bin 
(J MFdM/ J FdM). The systematic uncertainties are dominated by the uncertainties in 
the jet energy scale, which are 7% (30%) at 200 (1000) GeV/c 2 . The other uncertainties are 
due to the luminosity measurement (6.1%); luminosity matching at low mass (4.9%); the 
unsmearing correction, 0.5% (3%) at 200 (1000) GeV/c 2 ; the vertex cut (1%); and the jet 
selection cuts (1%). The total systematic uncertainty is given by the sum of the individual 
uncertainties in quadrature. The bin-to-bin correlations of the uncertainties are shown in 
Fig. H|0. 
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FIG. 1. d 3 <j/dMdr]idr]2 for | r/j et | < 1.0. The D0 data are shown by the solid circles, with error 
bars representing the ±lcr statistical and systematic uncertainties added in quadrature (in most 
cases smaller than the symbol). The histogram represents the jetrad prediction. 

The histogram in Fig. [I] is a prediction for the inclusive dijet mass spectrum from the 
next-to-leading (NLO) parton level event generator jetrad jHJ. The NLO calculation re- 
quires specification of the renormalization and factorization scales (/i = 0.5i?™ ax where i?™ ax 
is the maximum jet Et in the generated event), pdf (CTEQ3M |T4|), and parton cluster- 
ing algorithm. Two partons are combined if they are within TZ sep = 1.372., as motivated 
by the separation of jets in the data ||. Choosing an alternative pdf (CTEQ4M ||15|| . 
CTEQ4HJ fL5| , or MRS (A') [ll|) alters the prediction by as much as 25%, and varying \x 
in the range 0.25.E™ ax to 2i?™ ax alters the normalization by up to 30% with some M de- 
pendence. The CTEQ3M and MRS (A') pdf's are fits to collider and fixed target data sets 
published before 1994. CTEQ4M updates these fits using data published before 1996, and 
CTEQ4HJ adjusts the gluon distributions to fit the CDF inclusive jet cross section mea- 
surement ||. Figure |3] shows the ratio (Data — Theory)/Theory for the jetrad prediction 
using CTEQ3M with fi = 0.5E^ ax . Given the experimental and theoretical uncertainties, 
the predictions can be regarded as in good agreement with the data. The data are also in 
agreement within the given uncertainties with the cross section measured by CDF pf. 
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TABLE I. Dijet cross section for | r/j e t | < 1-0, and the ratio k(\ r/j e t | < 0.5)/k(0.5 < | r/j c t | < 1.0). 
High (low) systematic uncertainties are the sum in quadrature of the uncertainties from the ±1<t 
variations in the energy calibration, the unsmearing, the vertex corrections, luminosity matching, 
jet selection, and the uncertainty in the luminosity. 

Mass Bin (GeV/c 2 ) d 3 a / dM dr}\dr} 2 Ratio of Mass Spectra 

Bin Min. Max. Weighted ± Stat. Error Syst. Low Syst. High k(\ ry jct | < 0.5)//c(0.5 < | rj iet \ < 1.0) 
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FIG. 2. The correlations between systematic uncertainties in bins of dijet mass (see Table |) 
for 1 77jct | < 10. The correlations are calculated using the average systematic uncertainty. The 
discontinuities arise from the uncorrelated errors (adjacent to correlations of 1.0) and luminosity 
matching. 

In Table O, we show the x 2 resulting from a fit of theory to our data, using the full 
correlation matrix between different mass bins. The choice of pdf and renormalization scale 
is varied; all choices give reasonable probability. 



TABLE II. x 2 values calculated for various theoretical predictions for the dijet mass spectrum 
with |?7jet| < 10, and for the ratio of cross sections (15 degrees of freedom). 
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The ratio 


K( \fjjet | < 0.5)/«(0.5 < | 77jet 


< i.o), 


given in Fig. 


| and Table g 


, exploits the 



high correlation between uncertainties in the measurement of the dijet mass spectrum. The 
resulting cancelation of uncertainties leads to a systematic error of less than 8% for all 
M. The uncertainty in the theoretical prediction of this ratio is less than 3% due to the 
choice of pdf, and 6% from the choice of renormalization and factorization scale (excluding 
\x = 0.25E^p ax ). The x 2 values are shown in Table O. The predictions are in good agree- 
ment with the data, except for /i = 0.25i?™ ax which is excluded by the data (x 2 = 40.5, a 
probability of 0.04%). 
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FIG. 3. The difference between the data and the prediction (jetrad) divided by the prediction 
for | r/j e t | < 1.0. The solid circles represent the comparison to the calculation using CTEQ3M with 
[i = 0.5S™ ax . The shaded region represents the ±1<t systematic uncertainties. The effects of 
changing the renormalization scale and choosing a different pdf are also shown (each curve shows 
the difference between the alternative prediction and the standard prediction). 



The ratio of the mass spectra can be used to place limits on quark compositeness. A 
mass scale A characterizes both the strength of the quark-substructure coupling and the 
physical size of the composite state. Limits are set assuming that A ^> \f§ (where \f§ is the 
center of mass energy of the colliding partons), such that quarks appear to be point-like. 
Hence, the substructure coupling can be approximated by a four-Fermi contact interaction 
giving rise to an effective Lagrangian JTJ C = A(2ir / A 2 )(q L j fi q L )(q L j IM q L ), where A — ±1, 
and qi represents left-handed quarks. Limits are presented for the case where all quarks are 
composite, showing both constructive interference (A - for A = — 1) and destructive inter- 
ference (A + for A = +1). Currently there are no NLO compositeness calculations available; 
therefore, the pythia event generator is used to simulate the effect of compositeness. The 
ratio of these LO predictions with compositeness, to the LO with no compositeness, is used 
to scale the jetrad NLO prediction, as shown in Fig. |j. 

We employ a Bayesian technique |T7] to obtain from our data a limit on the scale of 



compositeness. Motivated by the form of the Lagrangian, a uniform prior is assumed in 
£ = 1/A 2 , and a Gaussian likelihood function P oc e~ x I 2 is used. The 95% confidence 
limit in A is determined by requiring that ftp(£')d? = 0.95. Since the ratio at NLO is 
sensitive to the choice of \i and pdf, each possible choice is treated as a different theory. The 
most conservative lower limits on the mass scale at the 95% confidence level are found to 
be A + > 2.7 TeV and A~ > 2.4 TeV. These limits are incompatible with the suggestion of 
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FIG. 4. The ratio of cross sections for | rjj et \ < 0.5 and 0.5 < | rjj et \ < 1.0 for data (solid circles) 
and theory (various lines). The error bars show the statistical and systematic uncertainties added 
in quadrature, and the crossbar shows the size of the statistical error. 



a compositeness scale A in the 1.5 to 1.8 TeV range found from earlier measurements || of 
the high Et jet inclusive cross section. 

In conclusion, we have measured the cross section for the inclusive dijet mass spectrum 
for If/jet | < 1-0 with M > 200 GeV/c 2 , and the ratio of cross sections for |^j e t| < 0.5 and 
0.5 < | ?7j e t | < 1.0, as a function of dijet mass. The data distributions are in good agreement 
with NLO QCD predictions. Models of quark compositeness with a contact interaction scale 
of less than 2.4 TeV are excluded at the 95% confidence level. 
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